An AC Stark Gradient Echo Memory in Cold Atoms 
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The burgeoning fields of quantum computing and quantum key distribution have created a demand 
for a quantum memory. The gradient echo memory scheme is a quantum memory candidate for light 
storage that can boast efficiencies approaching unity, as well as the flexibility to work with either 
two or three level atoms. The key to this scheme is the frequency gradient that is placed across the 
memory. Currently the three level implementation uses a Zeeman gradient and warm atoms. In this 
paper we model a new gradient creation mechanism - the ac Stark effect - to provide an improvement 
in the flexibility of gradient creation and field switching times. We propose this scheme in concert 
with a move to cold atoms (~ 1 mK). These temperatures would increase the storage times possible, 
and the small ensemble volumes would enable large ac Stark shifts with reasonable laser power. We 
find that memory bandwidths on the order of MHz can be produced with experimentally achievable 
laser powers and trapping volumes, with high precision in gradient creation and switching times on 
the order of nanoseconds possible. By looking at the different decoherence mechanisms present in 
this system we determine that coherence times on the order of 10s of milliseconds are possible, as 
are delay-bandwidth products of approximately 50 and efficiencies over 90%. 

PACS numbers: 42.50.Gy, 42.50.Hz 



I. INTRODUCTION 



The development of a quantum memory - a device 
that can store conjugate quantum variables - is cur- 
rently being driven by the field of quantum information 
processing, a field encompassing quantum computing 
and quantum cryptography. An ideal quantum memory 
for these applications would be 100% efficient, with long 
and controllable storage times, high storage bandwidths, 
and delay-bandwidth products, and faithful retrieval of 
the stored information. 

Their speed and lack of interaction with the environment 
make photons an ideal carrier for quantum information, 
unfortunately these same properties make them difficult 
to store. In recent years, however, much progress has 
been made towards the development of an optical quan- 
tum memory with techniques such as electromagnetically 
induced transparency (EIT), where recall efficiencies 
of over 40% in atomic ensembles jil] and storage times 
of over 1 second in solid state systems Q have been 
achieved. Another quantum memory candidate is the 
atomic frequency comb (AFC) scheme, based on photon 
echoes, where the re-alignment of atomic dipoles is 
required for the stored light to be re-emitted. Storage 
of coherent states with a mean photon number of 
approximately one has been demonstrated with AFC 
[3| and a maximum efficiency of 35% has been achieved 
for a fixed storage time on the order of microseconds 
(determined by the bandwidth of the system) 
The gradient echo memory (GEM) technique is also a 
photon echo based coherent memory. The key ingredient 
for GEM is a frequency gradient imposed along the 
memory, as illustrated in Figure [Ha). Not only does this 
gradient define the bandwidth of the system and cause 
re-emission of the pulse in the forwards direction by its 



reversal (Figure Hfb)), it also allows for 100% retrieval 
of the stored pulse if it is monotonic [5] . By altering the 
gradient, spectral manipulation of the pulse is possible 
[f|. In two level solid state systems an efficiency of 69% 
for GEM has been demonstrated Q. This work was 
carried out using an electric field (i.e. dc Stark) gradient 
created by placing four electrodes at the corners of the 
ensemble. 

With a two level system the storage time of GEM will 
be limited by the decay rate from the excited state. By 
moving to a three level A system in the far-detuned 
regime, as shown in Figure Uta), an effective two level 
atom is created. This is known as A- GEM, and apart 
from the increase in storage times achievable H, it 
also allows for pulse resequencing 0. The decay rate 
of this system is now determined by the ground state 
decoherence rate. 

Previous experimental work on A-GEM has used warm 
(65-70° C) gas cells, of length 7.5-20 cm and diameter 
2.5 cm, containing 87 Rb. This set-up has achieved 
efficiencies of up to 41% Q and coherence times of 20 /is. 
The frequency shift for the rubidium ensemble was 
created using a magnetic field (i.e. Zeeman) gradient, as 
alkali elements do not have a linear dc Stark shift. Using 
magnetic fields created by applying currents to solenoids 
wrapped around the gas cell, as in the above work, there 
is a lack of precision control over the gradient. Transient 
fields, which occur during gradient switching due to the 
inductance of the coils, also limit the switching time and 
can affect the rephasing process. 

One option for improving the gradient creation and 
control would be to move away from magnetic fields, and 
the coils necessary to create them, to an ac Stark (acS) 
shift. This would allow for an all-optically controlled 
quantum memory. Another option for improvement is 
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to move from warm to cold atoms. Due to the small 
decoherence rates in cold atomic ensembles this would 
allow for longer storage times and large on-resonance op- 
tical depths, due to the increase in density of the atoms. 
Implementing these improvements in concert would be 
beneficial as the acS effect is intensity dependent and 
cold atoms can be persuaded to occupy small volumes, 
reducing the laser power necessary. 

This paper investigates the feasibility of using an acS 
generated frequency gradient for A-GEM in an ensemble 
of cold 87 Rb atoms. After an overview of GEM theory 
in Section [ill the main body of the paper (Section IIIip 
discusses the proposed experimental implementation, 
including the theory behind the ac Stark shift and how 
it would apply to the creation of a frequency gradient 
across an ensemble of alkali atoms, as well as the optimal 
experimental parameters. Finally, factors that may limit 
this scheme will be discussed in Section IIVI such as 
coherence times and maximum efficiencies possible. 



II. GRADIENT ECHO MEMORY THEORY 
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FIG. 1: The GEM Scheme, (a) At time t = a pulse with 
envelope £(t) enters the atomic ensemble of two level atoms 
with ground state |1) and excited state \2) and with a linear 
frequency gradient T){z) = r\ x [z — L/2) applied across them. 
(b)(i) At time t — r the gradient is switched from 77 — > — r/ 
causing a rephasing of the atomic dipoles and (ii) a release of 
the pulse, time reversed, occurring at time t = 2r. 



Consider a collection of N two-level atoms with ground 
state |1) and excited state |2) with a resonant frequency 
uj and excited state decay rate 7, as shown in Figure 
IHa)(ii). We can then define atomic operators &ij = \i) (j\ 
and the atom-light coupling strength between the two 
levels 



where (i±2 is the dipole transition strength between the 
two levels due to the presence of a light field with an 
interaction volume V. 

If a linear frequency gradient is applied along an ensemble 
of atoms of length L, then the detuning from resonance of 
the atoms (i.e. the two-photon detuning §2p - see Figure 
2(a)(ii)) will be given by S2 P (z) = r/x (z — L/2). When a 
light field with a slowly varying envelope operator £(z,t) 
and centre frequency uj is sent into such an ensemble, 
the equations that govern the storage of the light, in a 
moving frame at the speed of light (t — ► t — z /c) and in 
the weak probe regime ((<5"ii) « 1), are [5j 

d t a 12 (z,t) = -[j/2 + ir](z-L/2)]& 12 (z,t) 

+ig£{z,t) (2) 



gN 

d z £(z,t) = i 0-12 CM). 

c 



(3) 



The efficiency of this writing stage of the memory is given 
by e w = 1 — exp(—2-Kd') [101 ] - determined by the effec- 
tive optical depth d! = g 2 N/(cg). This assumes that the 
bandwidth of the pulse is smaller than the bandwidth of 
the ensemble given by B s = r\L and also that B s ^> 7. 
To recall the pulse the gradient must be switched from 
1] — > —77. This causes a rephasing of the dipoles, and 
therefore a time-reversal of the initial storage process, re- 
sulting in the emission of a photon echo from the memory 
in the forward direction. If, as shown in Figure [TJ the in- 
put pulse enters the memory at time t = and the field is 
switched at t = r, then the output pulse will be released 
at t = 2t. Due to the reversal process the output pulse 
will be a mirror image of the input pulse with respect to 
time, i.e. £{z 7 t) — > £{z,—t) as shown in Figure [T{b)(ii). 
This can also be explained using the polariton descrip- 
tion of the storage process presented in (Tl| . 
If the gradient is monotonic along the memory then the 
read efficiency e r — e w to give a total read/ write effi- 
ciency for the memory of 



£rt« = [1 — exp(—2ird')] 



(4) 
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which will approach 100% for large optical depths 
d! —¥ 1. If the gradient is not monotonic the pulse 
will be partially re-absorbed as it leaves the memory, 
lowering the recall efficiency, with a maximum efficiency 
of 54% possible with no gradient [l2j]. Apart from the 
read-write efficiency there is also the storage efficiency 
e s to consider, which depends on the decay rate from 
the excited state and is of the form e s = exp(-jt). 

In previous work Q we have made use of the well 
known equivalence between a far detuned A system, with 
a strong coupling field and a weak probe field, and a 
two level system. This equivalence is illustrated in Fig- 
ure (Ha), showing the coupling field Rabi frequency fl c , 
probe field £ p , and the detuning for both detuned from 
the excited state Ai p (the one- photon detuning). The 
advantage of the A system (with g now dependent on 
/U13 for the probe field and [123 for the coupling field) 



3 




6.8 GHz 



FIG. 2: Atomic Level Structures. (a)(i) The three-level sys- 
tem used for A-GEM showing the decay rate 7 from the 
excited state |3), the decoherence rate 7<, between the two 
ground states |1) and |2), the one-photon detuning Ai p , the 
two-photon detuning 82? and the coupling strength between 
the two levels g. In the presence of a strong coupling field with 
Rabi frequency Q. c and weak probe field £ p this become equiv- 
alent to the two- level system shown in (ii), where the effective 
coupling strength g' = gQ c /Ai p . (b) The level structure of 
rubidium 87 showing the hyperfine splittings A^fs between F 
levels for both ground and excited states, as well as the fine- 
structure splitting Af s between the two excited states. Also 
shown are the detunings Aj/ F for a given laser frequency uji 
(see definitions in text). 



from resonance |Ai p | 3> dj; and (ii) Tjd ^> 1, where T is 
the fastest timescale of the system, usually either depen- 
dent on the pulse length t p or Q c . The equations of mo- 
tion for the A system then become equivalent to the two 
level equations of motion except for g — > g' = gCl c / Ai p . 
The effective optical depth for the equivalent two level 
system will therefore be 



d' 



g*N 

CI] 



Ai p 



(5) 



III. AC STARK SHIFT PROPOSAL 

The envisaged experiment is shown in Figure EJa). 
The probe and coupling fields would be created similarly 
to the method described in Q, with the coupling field 
being 6.8 GHz detuned from the probe to allow both 
S1/2 ground states of 87 Rb to be used (see Section IlII AD . 
There are three main components to this set-up: (i) 
cold atom storage; (ii) ac Stark gradient creation and 
(iii) switching (including probe and coupling field 
orientations). These will each be discussed in turn after 
firstly looking at the structure of rubidium, which will 
be used as the storage medium, and the ac Stark theory. 



A. Rubidium Structure 

Rubidium was used for previous experimental work as 
it provides a convenient working wavelength for the mem- 
ory. Rubidium 87, as with all alkali atoms, has a well 
known structure with, in this case, two hyperfine ground 
states Si/2,F = 1,2 and both a Dl and D2 transition, 
which correspond to P1/2 and P3/2 levels respectively 
with wavelengths of 795 and 780 nm [3] . This structure 
is shown in Figure [2jb) , containing the hyperfine struc- 
ture splittings Ahfs and A' h ^ s {J') for the ground and 
excited states respectively, as well as the fine structure 
splitting Ahf between the two excited states. 
The decay rate for the excited states J' of 87 Rb is j/2tt ps 
6 MHz, giving an excited state lifetime of approximately 
30 ns. Also, the hyperfine splitting Ahf s /2ir = 6.8 GHz 
for the ground states and A' hfs (3/2)/2n = 500 MHz, 
A' ft/s (l/2)/27r = 800 MHz for the excited states, with 
fine structure splitting Af s /2n = 7 THz. 



B. AC Stark Shift Theory 



over the two level one is that the decay rate from this 
excited state \2) is now limited by the decoherence rate 
7o *C 7- For large on-resonance, unbroadened, optical 
depths (d = g 2 NL/(ey) 3> 1) the conditions for this 
equivalence are [l3j : (i) the system being far-detuned 



When light of intensity I(r, z), frequency io\ and po- 
larisation q = 0, ±1 (corresponding to linear, right (+) 
and left (— ) circular polarizations respectively) is shone 
onto an atom there will be a change in energy of the 
internal states. This is known as the ac Stark effect. 
For alkali atoms, with structure as shown in Figure [^b), 
this effect can be calculated for a given ground state 
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FIG. 3: The Experiment, (a) Envisioned set-up for A-GEM 
experiment using cold atoms and an ac Stark gradient, with: 
Rb - atomic ensemble of rubidium 87; 50:50 - 50:50 beam 
splitter; PBS - polarizing beam splitter; PC - Pockels cell; 
BSh - Beam shaper; and HD - heterodyne detector. The po- 
larizations of the acS and trapping fields are assumed to be 
linear, while optimal probe and coupling field polarizations 
will be discussed in Section IIIIEI (b) Side-on view of the 
cylindrical atomic ensemble showing the direction of the in- 
tensity gradient I(z) and radius R. (c)(i) Gaussian and (ii) 
linear intensity profiles for the acS field per unit power over 
the ensemble. Here L has been taken to be 1 cm and R = 10 
/jm, determined by the radius of the trapping laser. 



\gi) = |l/2, F, m F ) from second order time-dependent 
perturbation theory to be [h| 



F,m,F 



2ce Q h ^ u>i — uj agi 

a. a L 



(6) 



where the sum is over all excited states \a) = \J',F', m' F ), 
uj agi is the frequency of the transition between |o) and \gi) 
allowing us to define the detuning as A = lo\ — Lo agi . This 
formula uses the rotating wave approximation, which is 
valid for detunings much smaller than the frequency of 
the transition, i.e. loi ~ io agi <C loi + & agi - For simplic- 
ity we can write the above equation as U Ftmp (A,q,I) = 
U FjmF (A, q)I(r, z), where U F . mF is the change in energy 
per unit intensity. An approximation to Equation [6j in 
the limit of A ^ A' h y s , is given by [It| 



U F , mF (A,q,I) 



irc 2 jl(r, z) 

2- 3 X 
2 + qg F m F + 1 - qg F m F \ 

^3/2, F ^1/2, F J 



where 7 and lo are the averaged values of the two excited 
levels, gF is the Lande factor [g\ = —0.5, 172 = —gi) and 
Aj^p is the detuning from the S 1 / 2 level (F = 1,2) to 
either J' = 1/2,3/2 levels, as illustrated in Figure^b). 
Apart from changing the energies of the atomic levels, 
there is the possibility that the atoms will absorb and 
then re-emit the light. This scattering of light by the 
atoms will affect not only the coherence time T co h achiev- 
able with the ensemble, but also the lifetime T trap of any 
trap that is used to contain them (see Section IIII CI) . 
The scattering rate for a given ground state will be de- 
termined by 



r F , mF (ui,q,i) 



I(r,z) 



E 



(9f\er ■ e qsc \a) (a\er ■ e q \g t ) 



LOl 



(8) 



again in the rotating wave approximation, where \gf) is 
the final state (which cannot be higher in energy than 
E gi + ftsLOi by conservation of energy) , w/, is the frequency 
of the transition from state \gi) to \gf) (negative if \gt) 
is lower in energy than \gi}), q sc is the polarisation of 
the scattered photon, and Lo agi is the frequency of the 
transition between states \gi) and \a). 
The above equation can be further simplified for alkali 
atoms if the detuning is much greater than the excited 
state hyperfine splitting A^ to give 



F,m f 



(ui,q,I) 



I(r,z) 



A, 



A 



LO fz ) 3 



3/2, fl4 



/2,F 



,(A,q)I(r,z) 



(9) 



where T F mF is the scattering rate per unit intensity, and 
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for all states \a) within the level J' . As can be seen from 
the above equations Tp mF « 1/A 2 while Up^ mF ^ 1/A 
and therefore, in the context of laser trapping, increasing 
the detuning for a constant trap depth will increase both 

T coh and Ttrap- 




1 10 100 1000 



-A/2jt (GHz) 

FIG. 4: Frequency Splitting, (a) Illustration of Sf splitting 
for the m,F = ±1,0 states of the F = 1, 2 levels. This is 
assuming a polarisation of q = 1 for the acS field. 5i t 2 is the 
difference between the splittings of the two mp =0 states of 
the respective levels from their original positions, (b) Log-log 
plot of (i) Si, (ii) &2 +Si, and (iii) St = &2 + 5i + <5i,2 as a 
function of detuning Ai/2, 2/271". The dashed line shows a fit 
to St for small detunings where 8 t oc 1/A. These traces were 
calculated using Equation [6] 



the F = 1 and 2 levels can similarly be found to be 

S h2 (A) = I (L\ (A) - c\ (A)) 

_ 7rc 2 7 / 2(A 3/21 - A 3/22 ) ^ A 1/2 ,i - A 1/2|2 \ 
2^3/i y A 3 / 2i2 A 3/2a A 1/2i2 A 1/2i1 / 

= ™*y A */' f ? + 1 W 

2w3/i \ A 3 / 2i2 A 3/2a A 1/2:2 A 1/2jl / 

These frequency splittings are illustrated in Figure HJJa) 
with the polarisation taken to be q — 1 (i.e. right circu- 
larly polarized) as, if linear polarisation were to be used, 
S F = 0. 

If, as was the case in 0, one rap state in each of the 
hyperffne levels is used for the ground states |1) and |2) 
from Figure [U(a), then we can define the total splitting 
in terms of Sf and <5i, 2 to be 

S t (A, q) = *i, a (A) - q (m 2 5 2 (A, |g|) + m^^A, |<j|)) . 

(14) 

Figure Hfb) shows the absolute value of 6t for m 2 = 
—mi = — 1 and the relative contributions from the three 
terms above as a function of detuning for q = 1. As 
can be seen, at large detunings Ai/ 2 ,f ^> Ahfs the split- 
ting between tuf levels becomes approximately equal (i.e. 
6\ ~ <5 2 ) and the relative contribution of <5i, 2 — > 0. This 
behaviour can be explained from Equations [T^] and 1131 
above as A rA ~ Aj/ :2 for A 1/2 ^ > A /l/s . 
From Equation [TJ] we can calculate the field gradient to 
be r)(z) ~ 2ir6 t d z I(r, z) assuming negligible change for 
the intensity in the transverse (r) direction. Using this 
equation for 77 allows the system bandwidth to be ex- 
pressed in terms of St and / as follows 

B s (A,\q\T) = ±-j\r,{A, q,I)\dz 

= \S t (A,q)(I(z = L)-I(z = 0))\.(15) 



Here we are interested in the energy splitting (and 
therefore frequency splitting hAv — AU) along an en- 
semble of atoms. This will determine the bandwidth of 
the system B s . The splitting per unit intensity between 
two adjacent tuf states in the same F level can be defined 
to be 



S F (A,\q\) = l/h\U FtQ (A,q)-U Ftl {A,q)\ 



(11) 



with the total magnitude of the splitting given by 
5f{A, \q\,I) = Sp(A, \q\)I(r, z). Combining the above 
equation with Equation [7] we obtain 
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3/2, F 



(12) 



The frequency splitting between the mp = states of 



C. Atom Trapping 

One of the most common methods of cooling atoms 
to millikelvin temperatures is a Magneto-Optical Trap 
(MOT) . It would seem an obvious suggestion to also use 
the MOT to confine the atoms during the memory pro- 
cess. However, a MOT has 6 circularly polarized beams 
detuned a few 7 below resonance and a magnetic field 
gradient on the order of 1 G/cm [l8l[T9j. The scattering 
rates associated with these beams and the mp splitting 
that would occur due to the magnetic field rule out the 
use of a MOT when storing information in an ensemble 
of atoms for A-GEM. If the MOT is turned off, however, 
the atoms are free to diffuse away from the interaction 
area and, though storage of light with a coherence time of 
1 ms has been achieved using this method [20| , if longer 
storage times are to be achieved another form of trapping 
must be used. 
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One way to achieve this trapping is to use a far off res- 
onance, red-detuned laser to create a dipole trap. This 
trapping field must be it polarized and detuned on the 
order of Af s or higher so as not to affect the different 
rciF splitting, and therefore the gradient, along or across 
the ensemble. Achieving tt polarisation requires the ad- 
dition of a constant dc magnetic field to create a quan- 
tisation axis. Such a beam, tightly focused and colli- 
mated, should provide sufficient trapping in the radial 
direction but much weaker trapping along the ensemble, 
as F t = -V(U t ) = -U t V(I(r,z)) and d r I > d z I. The 
radius of the ensemble will then become approximately 
equal to the waist of the trapping laser wt, which can 
be reduced down to the order of 10 /<m. This is desir- 
able for our system as a reduction in area for the acS 
gradient creating laser (A ac s > 2Lwt) will lead to an 
increase in the intensity, and therefore bandwidth, for a 
given power (see Section IIII D 21) . Realistic experimental 
parameters give L = 1 cm [2l|, |22j . To obtain the longest 
storage times possible the dipole laser should be detuned 
far from resonance. Assuming an easily obtainable ex- 
perimental wavelength of Xt = 1064 nm and a waist of 
Wt = 10 /mi, a trap of depth U t = 1 mK is achievable 
with a power of 1.5 W, giving a maximum scattering rate 

ofr t =4 S - x m. 

To increase the trapping in the z direction, one solution 
is to split the trapping field, using a 50:50 beam split- 
ter, and send it into the ensemble from both directions 
to create a standing wave trap, with trapping maxima 
occurring at At/2. The difference in detuning between 
two sites being 



2L 



(16) 



assuming a linear intensity profile. 

Dipole traps have already been used for pulse storage, 
with storage times up to a few milliseconds achieved [2J] . 
However, coherence times up to seconds are theoretically 
possible |25j |. 



D. Gradient Creation 

1. Wavelength Selection 

A critical parameter that must be determined is the 
wavelength of the acS laser to be used, as this will set a 
limit on the maximum frequency splitting possible for a 
given laser power and intensity distribution, as well as the 
scattering rate of the system. To optimize the wavelength 
we must balance the desired behaviour (i.e. frequency 
splitting) with the undesired effect of light scattering by 
the atoms determined by Tp iinF . It will be assumed that 
the acS field will be circularly polarized to maximize St 
for large detunings (see Section IIII Bp . 
From Equation we can see that the splitting, and 
therefore bandwidth, depends not only on 1/A but also 
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FIG. 5: The scattering rate required to produce a 1 MHz 
frequency splitting, for all m,F states of the F = 2 level, as a 
function of detuning Ax^^/^ (solid lines, left scale). Also 
shown is the intensity of the acS field required to produce 
the 1 MHz splitting (dashed line, right scale). These were 
calculated using Equations[6]and[8]without using the rotating 
wave approximation. 



on the ratio of A 1 / 2 .f/^3/2.f- As the detuning be- 
comes large compared to A/ s , A 1 / 2 .f/^3/2,f — > 1 as 
A3/2, f = A!/2,f + A/ s , and therefore the bracketed term 
in Equation 1121 heads to zero. This means that Sf will 
head towards zero faster than 1/A for large detunings. 
This can be seen from the difference in Figure [4jb) be- 
tween St and the dashed line, which shows the path St 
would take if it had only a 1/A dependence. This de- 
viation occurs on the order of THz detuning and is im- 
portant as it leads to Figure This shows the scatter- 
ing rates of different tuf states of the F = 2 level for 
a bandwidth of 1 MHz, assuming the same splitting ar- 
rangement used for Figure [U(b) and 5 = 1, i.e. 



2,m. 



■(A,l) 



10 6 f 2 , mF (A,l) 
**(A,1) : 



(17) 



with the corresponding intensity profile also shown [23j ]. 
From this figure it can be seen that the behaviour for 
the TOf = 0, — 1, — 2 levels is similar to that described for 
laser trapping in Section UlI Bl i.e. decreasing scattering 
rate with increasing detuning for a set trap depth. How- 
ever, for mp = 1 a minimum for the scattering rate is 
present at approximately 5 THz, with the scattering rate 
flattening out until 20 THz when it starts to decrease 
again. The mp = 2 state shows even more peculiar be- 
haviour, with a maximum appearing at approximately 
20 THz. A positive slope indicates that St is decreasing 
faster than F. This is to do in part with the dependence 
of St on A 1 /2 ) f/A 3 /2,f discussed above, as well as the 
relative strengths of A 3 / 2j gi and Ax/a )ffi and the levels 
which are allowed to contribute to T. All scattering rates 
converge at detunings much larger than Af s . 
The scattering rate for the F = 1 states are approxi- 
mately the same in magnitude but for the opposite mp 
state, with the mp = — 1 ground state containing this 
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minima. In both cases the unusual behaviour occurs 
for the states which are raised in energy with respect 
to mp =0. 

Due to the small probe approximation mentioned in Sec- 
tion |nl most of the population will remain in the F = 1 
state during the storage and retrieval processes. It would 
therefore be advantageous to make use of the minima 
and the F = l,m,F = —1 state as, though lower scat- 
tering rates can be achieved at much larger detunings, 
the laser intensities required become impractical. For in- 
stance, to achieve the same, or smaller, scattering rate 
than at the minima (approximately 11 s _1 ) requires a 
detuning greater than 40 THz. At this detuning, it re- 
quires nearly 10 times the laser intensity to achieve the 
same bandwidth as at the minima. For this minima de- 
tuning Si ~ <5 2 — 50 Hz/Wcm -2 and 5x,2 <C 5p. 



2. Beam Shaping 

The optimal intensity distribution - that which gives 
the maximum bandwidth for a given power P a - depends 
firstly on the orientation of the ensemble. As discussed in 
Section UlI C[ the ensemble will be cylindrical with L = 1 
cm and R — w t — 10 /mi. 

One of the simplest intensity profiles to produce is a fo- 
cused Gaussian, where 

I G (r',y)=I exp[-2r' 2 /wl], (18) 

with w the waist of the beam occurring at position y = 
(the centre of the ensemble, assuming the acS field is 
propagating along the y axis), I Q is the maximum inten- 
sity and r' is the radial component in the x-z plane. The 
only constraint on the acS beam shape is that it is mono- 
tonic along the ensemble (see Section HI|). In the case of a 
Gaussian beam, this means that over half the power will 
be lost, limiting the total splitting possible for a given 
laser power, with I a = 2P /(ttwI). 

To use the largest fraction of the beam possible, we can 
set w = 2L/3 so that if the centre of the beam is at 
2 = 0, approximately 99% of the remaining intensity will 
fall on the ensemble. Figure EJc) (i) shows this intensity 
profile along and across the ensemble for unit power. As 
the slope 77 is not linear, different frequency components 
will be stored with different efficiencies, as can be seen 
from Equation [5] In the case described above, where 
R <C L, the frequency change along the x direction for a 
set z position will be negligible. 

If one of the experimentally simplest intensity profiles for 
the acS beam is the one described above, then one of the 
most efficient will be one that covers only the ensemble, 
decreasing linearly from maximum intensity I at one end 
to at the other and with no change in intensity in the 
x direction. Such an intensity profile will be of the form 

I L (r',y) = I (1 - z/L) , \x\ < R 

= 0, |x| > R (19) 



over the length of the ensemble < z < L. In this case 
I = P /(LR) and therefore the maximum intensity, 
and bandwidth, achievable will be over 500 times larger 
than a Gaussian beam with the same power due to the 
smaller area the beam occupies. This can be seen from 
Figure G3c)(ii), which shows the intensity profile II for 
unit power along and across the ensemble. 
To change an initially Gaussian beam to one with an 
intensity profile like that from Equation [19] requires a 
beam shaper. These devices (for instance deformable 
mirrors, phase plates or liquid crystal spatial light 
modulators - LCSLMs) can be highly efficient (e > 0.9) 
and can be used to create nearly any desired beam shape 
with resolution on the order of 1000x1000 pixels for 
LCSLMs. This not only provides us with a method for 
optimizing the acS laser intensity profile, but would also 
allow for spectral manipulation of the pulse to be carried 
out with the ability to produce complex gradients and 
switching arrangements such as those described in [6j. 



E. Switching Protocols 

There are two components that make up the switching 
protocol: which mp states within the two F levels are 
to be used; and the method for switching the acS field. 
These will be discussed in turn. 



1. Selecting mp Levels 

Selection rules determine which mp states of the two 
different levels can be used, depending on the polariza- 
tions of the probe and coupling fields (q p and q c respec- 
tively), with the total change in angular momentum be- 
ing given by 

Amp — m% — mi — q p — q c . (20) 

As mi = —1 has already been decided upon as state 
|1) (see Section fill D 1|) . the above equation can be re- 
arranged to give 

m 2 = -1 + q p - q c . (21) 

We can determine the probe and coupling polarizations 
that will produce this maximum splitting by substituting 
the above relationship into Equation [14] to give 

St = Si,2 + 

q (<5i - (-1 + q P - g c )<5 2 ) 

~ qS F (2 + q c -q p ), (22) 

with the approximations S± ~ S2 ^> ^1,2 and \q\ 7^ 
discussed previously. This reveals that St can range from 
to 35p (remembering \m2\ < 2) with the latter being 
possible for either q c = 1, q p = or q c = 0, q p = — 1. 
Using the previous value of 8p determined in Section 
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MID 11 this would give a system bandwidth ol B s = 150 
Hz/Wcm~ 2 . If we combine this value with the maximum 
intensity achievable with the intensity profile II we find 
B s = 150 kHz/W and therefore system bandwidths on 
the order of 1 MHz would be obtainable with less than 
10 W of acS laser power using the optimal detuning. 
Also, using one of the above optimal level schemes with 
St = 3Sp, the scattering rate per 1 MHz splitting given 
in Equation [TTl which was calculated using St = 25 f, 
will be reduced as the intensity required to reach the 
same bandwidth will be reduced by a factor of 2/3. 
The polarizations for the probe and coupling fields mean 
that they cannot be completely separated on a polarizing 
beam splitter. To detect only the probe at the end of 
the memory a frequency selective measurement, such as 
heterodyne detection, could be used. 




FIG. 6: Switching Methods. (a)(i) Intensity profile per unit 
power and (ii) corresponding frequency splitting along the 
ensemble for switching method 1, with initial (solid line) and 
final (dashed line) gradients shown. Here the polarisation of 
the acS field q — 1 for both initial and final gradients. (b)(i) 
Intensity profile per unit power and (ii) corresponding fre- 
quency splitting for switching method 2, with initial (solid 
line) and final (dashed line) gradients shown. Here the initial 
polarisation of the acS field is q = 1, while the final polarisa- 
tion is q = — 1. This causes an frequency offset to the gradient 
of —St.max- Both (a) and (b) use the linear intensity profile 
II (see Equation ll9|) and the maximum St — 35 f- 



2. Field Switching Method 

To cause rephasing of the atomic dipoles we must be 
able to invert the detunings of the atoms. If we make 
the same assumption as in previous sections, i.e. |g| = 1, 
two switching methods become readily apparent. 
The first method involves reversing the intensity profile 
along the ensemble I(z) — > I(L — z). This is equivalent 
to the field switching method used in 0, [|], H| in that 



by reversing the intensity profile about the centre of the 
trap (z = L/2) the detunings are also reversed about 
this point, i.e. St (z) — > St(L — z). The intensity profiles 
for this method of switching are shown in Figure H)][a)(i), 
with the corresponding frequency gradients shown in 
Figure [Ha)(ii). This process involves no change in 
frequency of the stored pulse with respect to the input 
pulse. 

The second method involves switching the polarisation 
of the field q — > —q while keeping the same intensity 
gradient, as shown in Figure [5Jb) (i). This is a slightly 
more complicated process as the detunings are no 
longer reversed around the centre of the ensemble, with 
8t(z) — > — St(z). This method still results in an echo 
being produced. However, the stored pulse will now be 
frequency shifted with respect to the input pulse, as 
can be seen from the corresponding frequency gradients 
shown in Figure E^b) (ii) . This is because a switch from 
St(z) — > —5t(z) is equivalent to a switch about the centre 
and an offset of 5t, max being added. In a three level 
system this frequency shift can be overcome by altering 
the coupling field frequency in such a way as to cancel 
the initial shift. It should also be noted that Si. 2 will 
not contribute to St for this switching method as it only 
depends on the detuning and intensity of the acS laser, 
both of which will be constant. 

This second switching method would suggest itself as 
the easiest to implement as all that is required to switch 
rj would be a Pockels cell, which have switching times 
down to nanoseconds. It would not, however, allow 
for different frequency gradients and therefore filtering 
or manipulation of the pulse. Apart from this, if the 
second switching method is used, then the minima for 
the scattering rate found in Section MI D II will only 
apply for either the read or write stages of the memory. 
This is because, if the polarisation is reversed, the 
scattering behaviour of the levels will be reversed (i.e. 

f 2)1 ( g = -i)^r 2i _ x (g = i). 

The first switching method would allow for different 
gradients but involves much longer switching times (on 
the order of milliseconds for LCSLMs). The combination 
of beam simpers (BSh) and Pockels cells (PC) shown 
in Figure [3ja) allow for flexibility in beam shaping and 
fast switching times. If no spatial filtering is desired 
then only PC2 and BShl are needed, with the beam 
shaper determining the shape of the gradient and the 
polarisation switch causing the rephasing of the atoms. 
To allow different gradients to be used, an extra Pockels 
cell PCI and beam shaper BSh2 can be used. In this 
case, the second gradient can be prepared in advance 
and PCI used to select which beam shaper to use. The 
acousto-optic modulator (AOM) can be used to switch 
the acS beam on or off to decrease the scattering rate 
due to this field (see following section). 
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IV. LIMITING FACTORS 
A. Time Scales 

One main advantage of moving from warm to cold 
atoms is the extended storage times that can be achieved. 
There are two timescale of importance to our memory: 
(i) the trap lifetime T trap ; and (ii) the coherence time 

7~coh - 

The trap lifetime depends on both the scattering rate of 
the trapping laser, as well as the rate of inter-atomic col- 
lisions. These will also affect the coherence time. Using 
the trap parameters from Section [ill CI gives a coherence 
time of approximately 1/T t = 250 ms. Coupling the 
trap scattering rate with the recoil energy per emission 
E rec = (hk) 2 /2m PH, where m is the mass of rubidium 
and k is the wavevector of the transition (here taken to 
be the D2 transition), give the trap lifetime for a given 
depth U t to be 



7~trap 



mUt 
h 2 k 2 T t 



(23) 



Using this equation and trap parameters determined pre- 
viously (see Section |n| we find that the trap lifetime will 
be on the order of 10 seconds and should therefore not 
affect the coherence time of the system (see below) . 
Another effect that must be considered is the inelastic 
collision rate between the atoms in the trap, as this will 
also cause a loss of coherence as well as trap population 
if the atoms are not all situated in the F = 1 level. The 
collision rate in cold atoms has been studied in depth (see 
for example [H, I27W29I I ) with the collisional loss rate of 
the trap being given by 



dN 
~dt 



aN-p n 2 (r,t)d 3 r, 



(24) 



where n(r,t) is the density profile of atoms. As can 
be seen from the above equation there are two compo- 
nents to the decay: one due to background gas collisions 
with coefficient a being determined by the background 
gas pressure; and one due to inelastic collisions between 
atoms in the trap with coefficient j3. 1/a for a dipole 
trap can, in general, be approximated to 1 s for a trap 
pressure of 3 x 10- 9 mbar fijl. 

/3 is perhaps more interesting as it can give an indication 
of the time between inelastic collisions of the Rb atoms 
themselves. By using low trap laser intensities, the rate 
of hyperfine changing collisions has been estimated to be 
between fi hcc = lO" 11 ->■ 10" 10 em's" 1 [U H|. (3 hcc 
will be a lower bound on the total collision rate (as there 
are also non-hyperfine changing collisions to take into ac- 
count) but will be used to give an approximate rate of 
collisions. Here we will take j3h cc — 5 x 10~ n cm 3 s _1 . 
For the densities we are expecting (n sa 10 11 atoms/cm 3 ), 
and assuming a constant density within the trapping vol- 
ume for simplicity, Equation [23] will give an initial, and 
maximum, collision rate of j3n ~ 30 s _1 . 



Apart from the two sources of decoherence mentioned 
above, which will also affect the trap lifetime, there are 
two others that must be considered, namely the scatter- 
ing rate due to the acS field as well as the coupling field. 
The scattering rate for the acS field can be determined as 
a function of system bandwidth B s similarly to Equation 
[[7]to be 



r ac (A) 



l(A,1) 



MA) 



B 8 



(25) 



For the optimal detuning of Ax/2,2/^ = 5 THz discov- 
ered in Section UlID 1 1 for the F = I, nip = — I state with 
5=1, this simplifies to r ac ~ 7 x 10~ e B s for 5 t ~ 35f- 
The optical depth of the system depends on the Rabi 
frequency of the coupling field fl c and the one-photon 
detuning Ai p , as was shown in Equation [5] It therefore 
makes sense to express the effect of the coupling field on 
the ground state in terms of these two parameters by us- 
ing the relation between coupling field intensity I c and 
Rabi frequency I c = 2h 2 e cQ 2 / ii 2 3 to give 

r c (Ai p , n c ) = ^Lf2£f (A lp _ A hfSl q c ) n 2 c . (26) 

M23 

For pulses much longer than 1 fis and |fi c /Ai p | > 0.001, 
then T c > r QC . 

The total scattering rate for the system is simply the 
sum of the individual rates. The acS and coupling fields 
are only needed during the reading and writing phases 
of the memory process which each last for a minimum 
period determined by the pulse length, i.e. t r / w > t p , for 
single pulses but will become longer for multiple pulse 
storage. If we define a background scattering rate Tb g 
to include all the decoherence effects that are constantly 
present, i.e. scattering from the trapping laser, collisions 
and loss from the trap, and a read/write scattering rate of 
r rtu = r ac + T c , then we can determine the total storage 
efficiency to be 

e«(* P , T ) < exp (-2t p T rw ) exp (-(2t p + t s )T bg ) , (27) 

where t s is the time the pulse is stored in the memory. 



B. Efficiency 

There are three main experimental factors that will 
affect the storage efficiency of this system: the number 
of atoms N initially present in the MOT; the loading 
efficiency from the MOT to the dipole trap e^; and the 
size of the gradient which is applied \f](z)\ = 2ttB s /L. 
Substituting these values into Equations [5] and [4] gives 
the read/write efficiency of the system 



1 — exp — 




(28) 



To investigate the maximum efficiency as a function of 
storage time we assume, to maximize optical depth and 



10 




t u. 1000 
Ls/Lp 



FIG. 7: Memory Efficiency, (a) Total efficiency of the A- 
GEM memory as a function of pulse length t p for storage 
times of t a = t p and different ratios of |f2 c /Ai p | = (i) 0.01, 
(ii) 0.003, and (iii) 0.001. (b) Total efficiency for storage of 
pulses of length t p = 20 (is as a function of normalised storage 
time ts/tp for (i) single pulse; and (ii) multiple pulse storage, 
with |n c /Ai p | = 0.02. For all traces Ai„ = -2tt x 2 GHz, 
A ac = -2ir x 5 THz, B s = Bg = 9^2/(vrt p ), q c = 1, q p = 0, 
g = 2tt x 1.5 MHz, TV = 2.5 x 10 6 atoms, e L = 0.4, L = 1 cm 
and i? = 10 /^m. 



minimize scattering, that the bandwidth of the system is 
equal to the bandwidth of the pulse, i.e. B s = B p . For 
Gaussian pulses we can define the pulse bandwidth to 
be Bq = 9V2/ (nt p ), assuming 99% of the electric field 
is stored in the memory. Figure [TJ^a) shows the total ef- 
ficiency e t = e rw e s for different values of f2 c /Ai p , with 
storage time t s = t p , and using B s = Bq- 
As can be seen from the figure, for each ratio of f2 c /Ai p 
and storage time there is a maximum efficiency. This is 
a combination of two effects: firstly the increase in opti- 
cal depth as the pulse length increases due to the smaller 
bandwidths required, and extra decoherence that will oc- 
cur as the storage time increases. The decoherence is 
greater for larger values of f2 c /Ai p due to the increase in 
T rw this entails, and this can also be seen from Figure 
EJa). In this regime, however, the main source of de- 
coherence is the atomic collisions, with an approximate 



coherence time of 30 ms. 

For short pulses (t p -C 1 ms) and large values of 
|f2 c /Ai p | ~ 0.02 efficiencies approach unity as optical 
depths will be high and the background decoherence ef- 
fects will be negligible. This is therefore a good regime 
to investigate the delay-bandwidth product (DBP) of the 
memory. Here we define DBP= t s /t p as this will give an 
indication of the number of pulses (or bits) that can be 
stored in the memory at one time. 

Figure [7](b) shows e t for (i) single and (ii) multiple pulse 
storage with t p = 20 /its. The difference between the two 
is that, for multiple pulse storage, the coupling and acS 
fields are left on at all times, while this is not the case for 
single pulse storage. As can be seen, single pulse storage 
with efficiency greater than 90% (dashed line) occurs up 
to approximately t s = 130t p . For multiple pulse stor- 
age, this drops to 50 pulses due to the extra coupling 
field scattering, while 350 can be stored above the classi- 
cal efficiency limit of e t — 0.5. Currently a maximum of 
four pulse storage has been achieved with A-GEM with 
et — 1 [9]. To store pulses of 20 fis length requires a 
system bandwidth of approximately 200 kHz and there- 
fore an acS laser power of less than 2 W in the optimized 
regime. 



V. CONCLUSIONS 

We have made two proposals to improve the current 
experimental implementation of A-GEM. First, a new 
method of gradient creation using the ac Stark effect. 
Second, a move from warm to cold atoms, along with the 
longer coherence times this entails. 

This paper investigated the experimental viability of 
these proposals. Firstly, the trapping mechanism was 
investigated and it was determined that a far-detuned 
standing-wave dipole trap would not interfere with the 
acS gradient and have a lifetime on the order of seconds, 
while providing a small area for the acS laser. It was 
found that there is an optimal detuning for the acS field 
of approximately 5 THz to minimize scattering and max- 
imize bandwidth. With this detuning, and an optimized 
beam shape and level scheme, a bandwidth of 150 kHz/W 
could be created. Different methods for gradient switch- 
ing were investigated and a scheme devised that would 
allow a switch between any two arbitrary gradients on 
the order of nanoseconds using Pockels cells. 
Finally, factors that would limit the system such as scat- 
tering due to the trapping, coupling and acS lasers, as 
well as collisions between atoms were examined. These 
were combined with the read/ write efficiency to model 
the total system efficiency as a function of storage time 
and it was seen that for long pulses {t p 3> 1 ms), the 
coherence time was limited by the atomic collisions to be 
on the order of 10s of milliseconds. It was also found that 
for short pulses (t p <C 1 ms), efficiencies approach unity 
and 50 pulses could be stored in the memory at one time 
with efficiency greater than 90 %. The acS laser power 
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required to produce the necessary memory bandwidth to 
store these pulses would be less than 2 W in the opti- 
mized regime. 

We therefore conclude that using an ac Stark gradient 
with cold atoms is an experimentally viable option for im- 
proving the A-GEM system in terms of gradient creation, 
switching and manipulation, as well as storage times and 
time-bandwidth products achievable using cold atoms. 
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